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SUMMARY 

A  theory  was  recently  proposed  by  Wilson  and  Lin  to  explain  precursor  effects 

ahead  of  shock  waves.  An  indirect  test  of  this  theory  was  made.  Experiments  were 

performed  in  a  pressure  driven  snock  tube.  The  driven  gas  was  Argon  and  the 

shock  Mach  number  was  13.  2.  Small  amounts  of  H,  (2-75uHg)  were  added  to  the 
.  3  c 

Argon  (2  x  1C  uHg).  Precursor  profiles  were  recorded  as  a  function  of  the  H 

partial  pressure.  These  profiles  are  in  good  agreement  with  the  theoretical 
calculations. 
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SECTION  I 


INTRODUCTION 

In  a  recent  report  by  Wilson  and  Lin*,  a  theory  for  precursor  ionization  in 
shock  tubes  was  presented.  The  theory  was  founded  on  the  following  experimental 
results  obtained  by  Lederman  and  Wilson4-: 

(1)  the  precursor  is  the  result  of  photoioni'.ation,  the  radiation  coming  from  the 
shocked  gas,  and 

(<2)  the  precursor  profile  does  not  change  with  time  in  a  shock  fixed  coordinate 
system. 

As  discussed  in  Ref.  1,  a  one-step  photoionization  of  the  Argon  driven  gas  can¬ 
not  account  for  the  precursor.  The  radiation  for  frequencies  larger  than  the  first 
ionization  edge  frequency  of  Argon  is  effectively  absorbed  within  a  centimeter  ahead 
of  the  shock.  Two  alternative  explanations  have  been  suggested.  The  first  assumes 
a  multiple-step  photoionization  of  Argon;  resonance  radiation  first  excites  the  Argon 
driven  ga3,  and  continuum  radiation  completes  the  ionization  process.  The  second 
explanation,  the  one  recently  suggested  by  the  authors  of  this  report4*  and  others, 
is  that  the  precursor  is  due  to  photoionization  of  impurities  which  are  present  in 
the  Argon  driven  gas.  The  Argon  provides  the  high  temperatures  behind  the  shock, 
but  the  radiation  from  this  hot  gas  is  absorbed  by  impurity  atoms  in  the  driven  gas. 

This  second  explanation  is  favored  by  the  authors  of  this  report  for  the  following 
reasons:  In  the  experiments  of  Ref.  i,  the  Argcn  driven  pressure  was  changed  from 
1  mm.  to  10  mm.  without  substantially  altering  the  precursor  profile.  (The  shock 
Mach  numbers  did  vary  but,  since  then,  tests  have  been  performed  with  different 
driven  pressures,  keeping  the  Mach  number  constant,  with  essentially  tne  same 
result.  )  It  is  felt  that  if  a  multiple-step  process  in  Argon  were  responsible  for  the 
precursor,  a  ten-fold  increase  in  Argon  driven  pressure  should  at  least  have  pro¬ 
duced  a  substantially  sharper  drop  with  distance  ahead  of  the  shock  in  the  precursor 
profile. 
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With  these  experimental  ramifications  in  mind,  the  theory  of  Ref.  1  evolved. 

This  theory  can  be  applied  to  impurities  in  the  Argon  driven  gas.  It  is  necessary  to 
know  the  type  and  amount  of  each  impurity  and  absorption  cross-sectional  information 
for  these  impurities.  An  attempt  is  currently  being  made  to  calculate  the  precursor 
profile  due  to  residual  impurities  in  the  Argon  driven  gas. 

If  the  precursor  profile  were  found  to  be  sensitive  to  the  addition  of  small 
amounts  of  a  single  impurity,  an  indirect  test  of  the  theory  could  be  made.  This 
impurity  could  be  added  in  sufficient  quantities  so  as  to  "swamp  out"  the  contributions 
to  the  precursor  from  the  residual  impurities  in  the  driven  gas.  By  adding  varying 
amounts  of  this  impurity  to  the  driven  gas,  the  predictions  of  the  theory  could  be 
compared  with  experiment. 

The  precursor  profile  was  found  to  be  very  sensitive  to  the  addition  of  stroll 
amount  of  f .  -  lOuHg)  to  the  Argon  driven  gas.  Moreover,  the  absorption  cross- 
section  and  its  behavior  at  the  first  ionization  edge  of  H^,  is  accurately  known.  To 
test  the  validity  of  the  theory,  therefore,  experiments  were  performed  where  smail 
amounts  of  H ,  were  added  to  the  driven  gas.  The  predictions  of  the  theory  have 
already  been  set  forth  in  Ref.  1  for  an  H^,  impurity  (these  calculations  will  be 
modified  in  this  report)  and  comparisons  will  be  made  herein  with  experiment. 

The  authors  would  like  to  thank  Mr.  Tony  C.  Lin  for  his  help  with  the  computa¬ 
tions. 


SECTION  II 
THEORY 


The  precursor  profile  is  calculated  from  the  following  equation  : 

An  f  r°°  N(v)  Q-(v) 

n(x)  = - S_  - 1 -  exp  [  -  e  ]  d\» 

4-tU  *•  J  „  x 

+  ng  |  N(y)  Q^v)  Q(y)  Ei(  -»)dvj 


(1) 


where  p  =  Q(u)  x  and  Ei  is  the  exponential  integral.  The  symbols  in  this 
equation  represent  the  following: 

x  distance  ahead  of  the  shock 


n  (x)  electron  number  density 
A  cross-sectional  area  of  shock  tube 


n^  number  density  of  absorbing  impurity  molecules 
U  shock  speed 

Vi.  the  frequency  of  the  first  ionization  edge  of  the  impurity 
particle 

N(y))  the  number  of  photons  per  square  cm  per  second  per  unit  frequency  range 
emitted  at  frequency  v  at  x=0 
Q{(v)  the  cross-section  for  ionization 
Q  (v)  the  cross-section  for  absorption 


In  Ref.  1,  this  equation  was  applied  to  photoionization  of  H^.  The  justification 
for  using  a  black  body  source  for  N(v)  was  presented  in  this  report.  Also,  the  in¬ 
tegration  was  carried  out  from  V=v;  =  — 4  ev  to  '>  =  oo.  If  Argon  is  the 

H,  h 

driven  gas,  aa  is  the  case  for  the  present  investigation,  then  all  photons  with 

15.  7  e  ’ 

frequencies  higher  than  the  first  ionization  edge  of  Argon,  —  —  > 

will  be  absorbed  at  the  shock  front  and  will  not  penetrate  into  the  precursor  region. 


For  this  reason,  the  calculations  of  Ref.  1  were  repeated,  integrating  Eq.  (1)  from 
v  _  ^15^4 — ev_  to  v  =  ^ — —  .  Since  the  integration  interval  is  small,  average 

values  of  Q(v)  =  10  ^  cm^  and  Q.  (v)  =  4  x  10  cm^  were  used  in  this  interval.  ^ 


The  results  are  plotted  in  Fig  .  1.  Notice  that  the  same  behavior  in  the  precursor 


profile  with  varying  H^,  partial  pressure  is  still  maintained,  namely,  an  increase 
in  the  precursor  level  with  increasing  partial  pressure  and  an  eventual  decrease 
in  the  precursor  level  as  the  partial  pressure  is  further  increased. 
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In  Fig.  i,  a  comparison  is  made  between  the  precursor  profile  ao  calculated 
in  Ref.  1  and  as  calculated  by  the  method  outlined  above  at  an  partial  pressure 
of  50|jHg.  The  shapes  of  the  profile,  are  approximately  the  same,  although  the 
number  dsnsities  obtained  by  the  latter  method  are  lower  by  almost  an  order  of 
magnitude  than  those  obtained  by  the  former  method. 

"Peak  behavior"  for  the  present  calculations  (Fig.  1)  is  illustrated  ir  Fig.  3. 
The  number  densities  are  plotted  as  a  function  of  partial  pressure  at  a  distance 
of  x  =  50  cm  and  x  ■  100  cm  ahead  of  the  shock.  The  partial  pressures  at  which 
the  peaks  occur  as  a  function  of  x  will  later  be  compared  with  experiment  (Fig.  15). 

SECTION  III 

EXPERIMENTAL  PROCEDURE  AND  RESULTS 

A  pressure  driven  shock  tube,  described  in  a  previous  report^,  was  used  in 
this  investigation.  For  all  measurements  the  driver  pressure  was  800  psi  il^,  and 
the  driven  pressure  was  l  mm  Argon  to  which  prescribed  amounts  (from  Z-75|iHg. ) 
of  H^,  impurity  were  added.  A  continuous  record  of  the  shock  speed  was  obtained, 
using  a  microwave  technique  described  in  a  previous  report^.  In  spite  of  the 
addition  of  varying  amounts  of  impurity,  the  shock  Bpeed  at  the  test  station  was 
maintained  at  .4 i5  cm^sec  within  about  1%.  This  corresponds  to  a  temperature 
of  about  11,  000°K  behind  the  shock. 

The  precursor  profile  was:  obtained  using  a  cylindrical  probe,  .  875  inches  long 
and  .  040  inches  in  diameter,  aligned  with  the  flow  and  located  along  the  axis  of  the 
shock  tube  (Fig.  4a).  A  bias  of  -  6V  was  applied  to  the  probe  so  that  the  current 
to  the  probe  was  ionic  current,  and  this  current  was  determined  by  measuring  the 
voltage  drop  across  a  10K  resistor  connected  in  series  (Fig.  4b).  Typical  scope 
traces  are  shown  in  Fig.  5a  and  Fig.  5b.  The  low  sensitivity  beam.  Fig.  5a,  de¬ 
termines  the  time  of  shock  arrival  at  the  probe,  while  the  high  sensitivity  beam, 
Fig.  5b,  gives  the  precursor  distribution  ahead  of  the  shock. 
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TO  convert  ion  sa.nratio,,  current,  into  ion  number  d.n.itie.,  the  probe  wa. 

ral.bra.eu  -sing  a  microwave  flow-through  cavity^.  For  each  H,  parti,!  pree.ure 

a  etrgle  cavity  re.onanc.  was  recorded  (Pi,.  5c).  Thi.  determine,  the  number 

den.„y  for  one  di.tance  ahead  of  the  .hock.  The  number  den.ity  a.  a  function  of  x 

..  then  obtained  by  assuming  a  proportionality  between  ion  .aturation  current  and  ion 
number  density. 

The  raw  data  i.  plotted  in  Fig.  6.  The  ion  .aturation  current  i.  plotted  a.  a 
function  of  distance  ahead  of  the  shock.  The  amount  of  impurity  added  to  ,h. 
Argon  driven  ga.  appear,  a.  a  parameter.  The  precursor  densities  increase  with 
increa.ing  H,  partial  pre..ure  and  then  begin  to  decrease  a.  die  partial  pressure 
..  further  increa.ed  -  the  „me  behavior  that  i.  predicted  by  the  theory. 

The  curve,  in  Fig.  7  illustrate  this  "peak  behavior".  The  number  den.ities 
are  plotted  a.  a  function  of  H,  partial  pr.s.ure  for  two  fixed  distance,  ahead  of 
die  shock  ,50  and  100  cm,.  The  large  uncertainties  in  the  partial  pressure,  at 
which  the  peaks  occur  lead  to  Urge  "error  bar."  on  die  experimental  point,  i„ 

Fig.  15.  where  the  theory  is  compared  with  the  experiment. 

SECTION  IV 

COMPARISON  OF  THEORY  AND  EXPERIMENT 

In  Figs.  8  thru  14.  the  experimental  prolile,  are  compared  with  theory.  Each 
g  re  is  for  a  different  partial  pressure.  The  microwave  measurement  for  each 
H,  partial  pressure  is  indicated  by  a  circle.  At  all  pressures  the  shape  of  ,h.  pro- 
file,  predicted  by  the  theory  is  in  fairly  good  agreement  with  die  experimental 
profile  shape.  This  wUl  be  further  discussed  below  (Fig.  16). 

The  number  densities  given  by  theory  and  experiment  are  in  better  agreement 
at  the  higher  H,  partial  pressures.  I„  fact,  the  discrepancies  between  theory  and 
experiment  increase  as  th,  H,  partial  pressure  decrease,  (at  a  partial  pressure  o, 
fUHg.  the  experimental  number  densities  are  more  than  an  order  of  magnitude  higher 
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than  the  theoretical  values).  For  all  par:ial  pressures  the  experimental  values 
are  highe.  -ian  the  theory  predicts.  It  is  felt  that  the  reason  for  the  higv  *r  experi¬ 
mental  values  is  that  the  assumption  of  no  reflection  of  the  radiation  by  the  shock 
tube  walls  *  is  not  entirely  valid.  Even  at  these  high  U.V.  frequencies  enough 
radiation  is  reflected  baoc  into  the  shock  tube  to  explain  the  discrepancies.  This 

has  been  experimentally  verified  by  measuring  the  precursor  ahead  of  the  shock  in 

4 

a  large  dump  tank  where  the  walls  have  been  removed. 

There  is  no  clear-cut  explanation  for  the  larger  discrepancies  at  the  lower 
partial  pressures.  It  is  very  likely  connected  with  the  fact  that  the  contributions  to 
the  precursor  from  other  impurities  are  more  evident  when  leBS  is  added  to  the 
driven  gas.  At  least  this  explanation  is  consistent  with  larger  experimental  than 
theoretical  values  of  number  density. 

In  Fig.  13,  (H^  partial  pressure  is  50pHg)  the  precursor  profile  obtained  from 
5 

Wetzel's  theory  ,  applied  to  photoionization  of  H^,,  is  plotted  'or  comparison.  The 
number  densities  for  tli'  Vetzel  curve  are  several  orders  of  magnitude  higher  than 
the  experimental  measurements,  and  the  "fall-off"  with  distance  ahead  of  the  shock 
is  not  as  fast  as  the  experiment  indicates. 

In  Fig.  15,  a  comparison  is  made  of  the  theoretical  and  experim  mtal  values  of 
partial  pressure  at  which  thj  peak  value  in  the  precursor  number  density  occurs 
(as  a  function  of  distance  ahead  of  the  shock).  These  values  of  partial  pressure  are 
obtained  from  curves  such  as  those  plotted  in  Figs.  3  and  7.  The  agreement  between 
theory  and  experiment  is  about  as  good  as  can  be  expected. 

The  precursor  profiles,  both  experimental  and  theoretical,  exhibit  exponential 
decay  far  from  the  shock  front.  In  Fig.  16>  the  slopes  at  a  distance  of  x  =  100  cm. 
are  compared.  The  slopes  are  obtained  directly  from  Figs.  8  thru  14.  The  theory 
accurately  predict  the  "profile  shape"  for  all  H^,  partial  pressures. 
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It  should  be  mentioned  that  such  factors  as  black  body  efficiency,  reflection 
of  radiation  from  the  walls  of  the  shock  tube,  values  one  chooses  for  Q(M),  Q.  (v), 
etc.  ,  affect  the  absolute  values  one  obtains  for  precursor  number  dcnsitiec,  but  the 
shape  of  the  profile  itself  is  not  very  sensitive  to  these  factors.  A  more  important 
criterion  for  the  efficacy  of  the  theory  is,  therefore,  how  well  it  predicts  the  profile 
shape  and  general  qualitative  features  of  the  profile  behavior  with  changing  impurity 
concentrations,  rather  than  how  well  it  predicts  the  absolute  number  density  levels, 
provided,  of  oo:  se,  that  the  values  are  not  too  far  off. 

SECTION  V 

CONCLUDING  REMARr  S 

The  main  objective  of  the  experiments  was  to  test  the  validity  of  the  impurity 
photoionization  theory  proposed  in  Ref.  1.  The  good  qualitative  and  quantitative 
agreement  between  the  theory  and  the  experiments  justiii »r  increased  confidence  in 
the  use  of  this  theory  for  calculating  precursor  profiles  under  different  operating 
conditions. 

One  such  application  might  very  well  be  to  determine  the  precursor  distribution 
ahead  of  a  reentry  vehicle  as  a  function  of  height,  vehicle  attitude  and  velocity,  etc. 
Since  the  composition  of  the  atmosphere  as  a  function  of  height  is  fairly  well 
established,  it  should  be  possible  to  make  tl  is  type  of  calculation. 

Further  experiments  are  planned  to  determine  the  types  of  ions  and  the  number 
densities  of  each  type  in  the  precursor  region.  This  kind  of  information  would 
prove  invaluable  in  tying  together  all  aspects  of  the  precursor  problem. 
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FIG.  4  ELECTROSTATIC  PROBE  AND  CIRCUIT 
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Typical  Experimental  Traces 

(a)  Low  Sensitivity  Electrostatic  Probe  Trace 
BIAS  =  -6V,  2V/DIV. 

(b)  High  Sensitivity  Electrostatic  Prcbe  Trace 
BIAS  =  -6V,  50  MV/OIV. 


(c)  S-Band  Cavity  Kesponse 
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FIG.  7  PRECURSOR  NUMBER  DENSITY  VS.  Hz 
PARTIAL  PRESSURE  -  EXPERIMENT 
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